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Abstract. - We study the phase dynamics in coupled Josephson junctions describing by system 
of nonlinear differential equations. Results of detailed numerical simulations of charge creation 
in the superconducting layers and the longitudinal plasma wave (LPW) nucleation are presented. 
We demonstrate the different time stages in the development of the LPW and present results 
of FFT analysis at different values of bias current. The correspondence between the breakpoint 
position on the outermost branch of current voltage characteristics (CVC) and the growing region 
in time dependence of the electric charge in the superconducting layer is established. The effects 
of noise in the bias current and the external microwave radiation on the charge dynamics of the 
coupled Josephson junctions are found. These effects introduce a way to regulate the process of 
LPW nucleation in the stack of IJJ. 
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Introduction. — The intrinsic Josephson junctions 
(IJJ) in high temperature superconductors (HTSC) are 
formed from the atomic scale superconducting layers (S- 
Iciyers), with each layer having a thickness that is com- 
parable with the Debye screening length. [1] Because of 
it, the electric charge docs not screen in S-layers perfectly 
and it leads to the capacitive coupling between junctions 
and the existence of the LPW along the stack. 
■ The inductive coupling [S] of IJJ in the absence of mag- 
iietic field can be neglected for the small size stacks, and 
the phase dynamics in such nanojunctions is determined 
by the capacitive coupling only. In this case the sys- 
tem of IJJ is described by capacitively coupled Josephson 
junctions (CCJJ) model or the model with the diffu- 
sion current (CCJJ-I-DC) [6j[7|- The investigation of the 
charge dynamics in such stack allowed us to predict new 
physical properties of the coupled system of JJ. IHIH] Its 
CVC is characterized by the multiple branch structure and 
branches have a breakpoint and some breakpoint region 
(BPR) before the transition to another branch. [3] 

The physical properties of IJJ arc investigating very in- 
tensively today. The collective Josephson plasma reso- 
nance was observed experimentally some time ago in Ref. 
in] but the observation of the powerful coherent radia- 
tion from the stack of IJJ |10| and the experimental man- 



ifestation of the breakpoint and the BPR [TT] in CVC 
of Bi2Sr2CaCu20y stimulates new investigations in this 
field. The key problems are the understanding of the 
mechanism of this radiation and the way to increase its 
power. An important question concerns the nucleation of 
the LPW in the stack of IJJ which is not investigated yet. 
The correspondence between the breakpoint's position in 
CVC and the parametric resonance region in time depen- 
dence of the electric charge in the S-layers is not precisely 
established as well. The important open issue is a question 
if it is possible to regulate the process of LPW nucleation. 

In this paper we investigate the phase dynamics in the 
coupled system of Josephson junctions. We show the dif- 
ferent time stages of charge creation in S-layers and the 
development of the LPW. The data concern the charge 
distribution along the stack and the results of fast Fourier 
transformation analysis at different values of bias current. 
It is found that the nucleation process of LPW can be 
affected by the noise in bias current and the external mi- 
crowave radiation. It makes possible to regulate the pro- 
cess of LPW nucleation in HTSC. 

CVC and time dependence of the charge in su- 
perconducting layer. — To simulate the CVC of IJJ, 

we solve a system of dynamical equations d^Lpi/dr"^ = 
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{I -siinpi - Pdipi/dr) + a{smipi+i +sin(^/_i — 2sinc^i) + 
af3{dipi+i/dT + dipi-i/dr — 2dipi/dT) for the gauge- 
invariant phase differences ipiir) = 9i^i{t) — 0i{t) — 
^ Ji~^^ dzAz{z,T) between superconducting layers (5- 
layers), where 9i is the phase of the order parameter in 
the S- layer /, is the vector potential in the barrier and 
a is the coupling parameter. The third term in the right 
hand side is the diffusion current which is determined by 
the difference of the generalized scalar potentials of super- 
conducting layers. [2lll2j In our simulations we measure the 
voltage in units of Vb = hujp/{2e) and the current in units 
of the Ic- We use a dimensionless time r = tujp, where ujp 
is the plasma frequency ujp = y^eic/^Z?, Ic is the critical 
current and C is the capacitance. The CVC and time de- 
pendence of the charge oscillation in the S-layers are simu- 
lated at coupling parameter a = 1, dissipation parameter 
/3 = 0.2 and periodic BC. We chose arbitrary initial value 
of bias current to record the time dependence. The time 
dependence consists of time and bias current variations. 
At each current step the phase dynamics is developed in 
the definite time interval. A small noise with its maximum 
in the interval (~10~*, -|-10~^) is added in bias current in 
our simulations. The details concerning the model and 
numerical procedure are presented in Refs. [HE] 




Fig. 1: (Color online) The enlarged part of the outermost 
branch of the CVC of the stack wit 9 IJJ. The inset shows 
the total CVC. 

To compare our results with previous investigations of 
the CVC and charge dynamics in the coupled system of 
JJ we concentrate here on the stack with 9 JJ. A corre- 
spondence between the CVC and charge oscillations in the 
S-layers in such stack was demonstrated in Ref. [S] It was 
shown that at chosen values of a and /3 the LPW with 
wave vector = Stt/Q is created in the stacks with 9 cou- 
pled JJ at the breakpoint. [4] Here we study the nucleation 
process of LPW which was not touched before. 

The CVC of the stack with 9 coupled JJ is presented in 
the inset to Fig. [TJ It has a large hysteretic region with 9 
branches. We show by circle the BPR on the outermost 
branch and its enlarged part in this figure. The arrows 



indicate the bias current values at which the detailed study 
of the charged dynamics have been done. We investigate 
the dependence of the electric charge in S-layer vs time at 
the following values of bias current / = 0.59, 0.572, 0.5714 
and / — 0.5711. We have done the FFT analysis of the 
charges oscillations at these values of current as well. As 
we mention in Ref. [5], the positions of the breakpoint in 
CVC and the onset of charge growth in S-layer in time 
dependence are not coincide. 

Time stages in the nucleation of LPW. A cor- 
respondence between the breakpoint in CVC and charge 
dynamics in S-layers was demonstrated in the previous 
research. [8] Here we demonstrate that the creation and 
development of LPW has different stages in time (or cur- 
rent). We show that the breakpoint in CVC reflects just 
the region of sharp increase of the LPW amplitude. Inves- 
tigation of the time dependence of electric charge in the 
S-layers give us the information concerning the nucleation 
of LPW. 

At chosen values of a and /?, the LPW with wave vector 
fc = 87r/9 is created at the breakpoint. [U Below we show 
the profiles of the charge oscillations and some enlarged 
parts for different values of bias current indicated in Fig.[T] 

Fluctuation region. Fig. [5^1 shows the charge-time de- 
pendence for the first two S-laycrs of the stack with 9 
coupled JJ at / = 0.59, which is far enough from the 
breakpoint in the outermost branch of CVC. We see that 
the charge value Qi/Qo is practically the same in both lay- 
ers and it is about 10^®, i.e. it is on the noise level. The 
charge oscillations here have mostly an irregular charac- 
ter. But nevertheless, as we can see in Fig. [2Jd, there are 
the regular oscillations in a short time interval with the 
frequency corresponding to the LPW. We call this part in 
CVC and the corresponding part in time dependence of 
the charge on the S-laycrs as a fluctuation region. 




Fig. 2: (Color online) (a) The characteristic part of dependence 
the charge in first two S-layers vs time at / = 0.59 for the stack 
with 9 coupled JJ; (b) The enlarged part of the time interval 
(140,190) which demonstrates the fluctuation region. The thick 
line shows the charge in the first layer of the stack, the thin 
line corresponds to the charge in the second layer. 

Let us clarify the shape of the charge signal we sec in 
Fig. [5)d . To do it, we demonstrate in Fig. [3] the time 
dependence of the voltages in the first Vi and last Vn 
junctions (left y axis), which are the adjacent junctions 
to the first layer (we use the periodic BC and first layer 



p-2 



Title 



has number 0). The charge density Qi in the S-layer I 
is expressed by the voltages Vi and V/_i in the adjusting 
junctions Qi = Qoa{Vi - Vi-i), where Qq = ££oVo/r\,, 
and is Debye screening length. On the right y axes 
wc plot the value of charge in the first S-layer, which is 
actually cc^ual to Vi — Vjq at a = 1. At some time in 
during the oscillations the sign of charge in the first layer 
is changed: the vertical dashed lines show the sign changes 
of the difference Vi ~ Vn . 
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Fig. 3: (Color online) The formation of the charge signal's 
shape in the first S-layer of the stack with nine IJJ. 



We see that at time moment corresponding to the first 
dashed line, the curve Vi is getting above of curve Vat, but 
at time moment corresponding to the second dashed line 
the situation is opposite and the charge in the first layer 
is getting negative. The time dependence of the voltage 
difference, i.e. charge in the first S-layer is shown by thick 
line. We see that the shape of the charge signal coincide 
qualitatively with the charge signal observed in Fig. ^jp. 
Of course, the shape of charge signal might have some de- 
viations from this one in the different layers. Its change 
in time reflects the character of Josephson junctions phase 
dynamics at different values of bias current. So, by shift- 
ing of voltages along the time axis in the Fig. [3l we may 
reproduce different shape of the charge signal ( particu- 
larly that one which appears in Fig. ^jp). 

The spectrum of Fast Fourier Transformation (FFT) 
analysis of the charge oscillations in the first S-layer in the 
fluctuation region is shown in Fig. As it's expected, it 
is broad and degraded. But nevertheless it demonstrates 
the manifestation of the LEW (maximum in the inter- 
val 0.22-0.25) and the Josephson oscillations (maximum in 
the interval 0.43-0.48). Of course, the charge distribution 
along the stack deviates from the wave with k = Stt/Q, as 
we can see in Fig.HJj. The solid line shows the sin(87r/9)z, 
the rhombus show the charge value in the flrst S-layer 
Qi/Qa at some flxed time moment t = 170.1. 
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Fig. 4: (Color online) (a) Results of FFT analysis of the charge 
oscillations in the stack with nine IJJ at bias current / = 0.59; 
(b) Charge distribution along the stack at / = 0.59. The sohd 
line in (b) shows the wave sin{8n 



Island region. With approaching the breakpoint on 
the outermost branch, the number of such intervals of reg- 
ular oscillations we demonstrated in Fig. ^jp is increasing 
and there size along time axes as well. The most important 
fact in this case is the value of the oscillations amplitude. 
It exceeds the noise value 10~^ essentially now. So, we call 
such regions like the region shown in Fig. [5^ at / = 0.572 
as the " island regions" . 

Fig. [S^ shows the charge-time dependence for flrst two 
S-layers of the stack in the island region and Fig. [SJd shows 
the enlarged part marked by circle in this region. We see 
that the amplitude of charge oscillations in shown island 
might have a value up to 7 * 10^^. 
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Fig. 5: (Color online) (a) The characteristic parts of charge 
time dependence in the stack with 9 IJJ in island region; (b) 
The enlarged part of time dependence marked by circle in figure 
(a). The thick line shows the charge in the first S-layer of the 
stack, the thin line corresponds to the charge in the second 
S-layer. 

In the island region we can precisely determine the LPW 
frequency by the FFT analysis in contrast to the fluctua- 
tion region. In Fig. |6ji we demonstrate the results of FFT 
analysis for charge oscillations at / = 0.572. As we can 
see, there is still the broadening part in the FFT spec- 
trum, but the charge distribution along the stack follows 
the modulated LPW with A: = 87r/9 (see Fig.Eb ). 

Alternative amplitude region. With increase in time 
(decrease in current) the islands are growing, joining, 
forming the region of charge oscillations with alternat- 
ing amplitude. The characteristic part of it is shown in 
Fig. [7^. The amplitude of oscillation here is around 10"'', 
i.e. one order of value exceeds the noise level. 
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Fig. 6: (Color online) (a) Results of FFT analysis of the charge 
oscillations in the stack with nine coupled JJ at bias current 
I = 0.572. (b) Charge distribution along the stack at / = 
0.572. The solid line in (b) shows the wave sm(87r/9). 
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Fig. 7: (Color online) (a) The characteristic parts of charge- 
time dependence in the stack with 9 IJJ in the alternative 
amplitude region; (b) Results of FFT analysis of the charge 
oscillations at bias current / = 0.5714. 
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Fig. 8: (Color online) (a) The characteristic parts of charge- 
time dependence in the stack with 9 IJJ in growing region; (b) 
The enlarged part of (a) in time interval 189800-189850. The 
thick line shows the charge in the first S-layer of the stack, 
the thin line corresponds to the charge in the second S-layer; 
(c ) Results of FFT analysis of the charge oscillations at bias 
current / = 0.5711. 



time development of LPW: the fluctuation region (LPW 
in short time interval on the noise level); the island re- 
gion (with charge oscillations exceeding the noise level); 
the alternating amplitude region (all islands are joined), 
and the growing region (charge oscillations with growing 
amplitude) . 



In Fig. [7}d we show the enlarged part of charge oscil- 
lations in the time interval (186900, 187000). We can 
see clear here that the charge amplitude in the second 
layer is smaller than in the first one. It shows that the 
LPW is not the 7r-mode, where k = n/d and the ampli- 
tude is the same in all layers. As we mentioned above, 
the LPW with k = 87r/9(i is created in the stack with 9 
coupled JJ at a = 1 and /3 = 0.2. In Fig. [T]; we demon- 
strate the results of FFT analysis for charge oscillations 
in the first S-layer at / = 0.5714. The LPW frequency is 
wlpvf/SttWp = 0.22583 which is practically twice smaller 
than the Josephson frequency ujj/2'nujp ~ 0.45166. 

Growing region. In Fig. [8^ and Fig. [8Jd we demon- 
strate the charge oscillations and some enlarged part 
for the beginning of the region with growing amplitude. 
Fig. [5]3 demonstrates clear the fact that the amplitudes 
of the charge oscillations in the neighboring layers in the 
growing region is essentially different. Fig. [St shows the 
results of FFT analysis of the charge oscillations in the 
first S-layer at bias current / = 0.5711. The peaks are 
very narrow and sharp as it should be at the resonance 
condition. In this region the Josephson frequency is equal 
to the LPW frequency very precisely. 

So, we may distinguish clearly four different stages in 
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Fig. 9: (Color online) Time dependence of the charge oscilla- 
tions in the first S-layer of stack with nine IJJ at different values 
of the noise amplitude in bias current. The insets (a) and (b) 
show the enlarged beginning parts of the growing region. 

Effects of noise level and external radiation on 
the LPW nucleation. — Let us briefly discuss the pos- 
sibility of affecting the process of LPW's nucleation. First 
we present the results concerning the influence of the level 
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Fig. 10: (Color online) (a) Effect of radiation on the LPW 
nucleation in the stack with 9 IJJ at A = 0; (b) The same at 
A = 0.3. 



of noise in bias current. As it was mentioned above, we 
add a small noise with amplitude 10~^ in bias current in 
our calculations. The noise stimulates the appearance of 
difference of the phase differences in neighbor junctions. 
It leads to some distribution of the charge in the S-layers 
along the stack. At current values much larger than I{,p 
this distribution of charge attenuates in time. Approach- 
ing the breakpoint, the corresponding oscillations are in- 
creasing, especially in the growing region. 

Fig. ini demonstrates the time dependence of the charge 
in the first superconducting layer of the stack with nine 
coupled JJ in the beginning of the growing region for three 
values of the noise amplitude: 10~*, 10~^° and 10"^'^. It 
shows a remarkable fact that the noise with larger ampli- 
tude creates the LPW at smaller time (current) values, 
than the noise with the smaller amplitude. The inset (a) 
allows us to see clear that the growing region starts first 
in case of the noise amplitude equal to 10~* and the inset 
(b) demonstrates that for the noise amplitude 10^^'^ the 
growing region starts in the last turn. This result is un- 
derstandable: it's easier to arrange a difference between 
neighbor junctions by larger value of the noise amplitude. 

The another possibility to affect the nucleation process 
is the effect of the external radiation. If we irradiate a 
stack of Josephson junctions with microwave of frequency 
LVr, it causes an alternating current through the stack. 
This alternating current can be described as the additional 
applied current Asm{LOrt), where A is a dimensionless am- 
plitude. 

In Fig. [TU] we demonstrate the effect of radiation on the 
charge oscillations in the first S-layer in the beginning of 
the growing region. We see that the radiation with fre- 
quency uj/uip = 0.5 and amplitude A = 0.3 (Fig. [TOb) 
shifts the starting point of growing region (shown by ar- 
row) in comparison with case A = (Fig.fTUk). Our results 
show that the change in the frequency leads to the shift of 
this point and additional resonance features. This effect 



of radiation will be considered in detail somewhere else. 

We would like to mention here that the nucleation pro- 
cess can be affected by change in temperature as well. 
The change in the McCumber parameter should lead to 
the creation of the LPW with another wave number [3] 
and, correspondingly, to the nucleation process specific to 
this wave number . 

In summary, we studied the nucleation of the longitu- 
dinal plasma wave in the coupled system of IJJ in HTSC 
and found the different stages in its process of develop- 
ment. The answer to the question concerning the corre- 
spondence between the breakpoint's position in CVC and 
the parametric resonance region in time dependence of the 
charge on the S-layers is found. We showed that the posi- 
tion of the breakpoint on the outermost branch of CVC is 
related to the region with sharp increase of the amplitude 
of charge oscillation in the superconducting layers. We 
demonstrated that the onset of the growing region can be 
shifted by noise in the bias current and microwave radia- 
tion. These effects open a way to regulate the process of 
LPW nucleation in the stack of IJJ. 
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